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I.  INTRODUCTION 


The  shielded  suspended  substrate  line  (SSSL)  is  a  transmission  medium  useful 
for  radar  and  circuits  at  microwave  and  millimeter-wave  frequencies.  In  order  to 
use  this  transmission  medium  in  the  construction  of  microwave  circuits  and  filters, 
it  is  necessary  to  have  valid  circuit  models  for  typical  discontinuities  such  as  the 
series  gap  in  line,  open-ended  stub,  and  a  discontinuous  change  in  width. 

There  is  a  definite  need  for  an  accurate  full-wave  analysis  of  strip 
transmission  line  structures.  Even  though  many  studies  have  been  performed  for 
the  discontinuities  of  conventional  microstrip  line,  almost  no  discontinuicy 
calculations  are  available  for  SSSL.  By  full-wave  analysis  is  implied  the  process 
of  rigorously  solving  the  frequency-dependent  electromagnetic  (EM)  boundaiy' 
value  problem  with  retention  of  all  the  field  components. 

In  this  thesis,  the  discontinuity  structure  of  interest  incorporated  in  a  SSSL 
resonator  has  been  calculated  by  a  fomiulation  employing  full-wave  analysis.  The 
solution  of  the  problem  has  been  derived  using  an  efficient  method.  Specifically, 
the  derivation  of  the  characteristic  equation  for  resonant  frequencies  of  the 
resonator  model  is  carried  out  using  Galerkin’s  technique  applied  in  the  spectral  or 
Fourier  transform  domain  instead  of  the  space  domain.  The  resonant  frequencies 
of  the  structures  of  interest  are  obtained  by  numerically  solving  the  characteristic 
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equation.  The  details  of  the  analysis  method  will  appear  in  Section  II  of  this 
work. 

In  Section  III,  the  step  discontinuity  will  be  analyzed  by  placing  the 
discontinuity  in  an  open-ended  strip  line  resonator.  In  the  continuous  strip  case,  the 
strip  line  has  no  discontinuity  step  in  width.  Following  this,  the  equivalent  circuit 
of  the  discontinuity  will  be  deduced  from  the  numerical  resonance  data  for  the 
embedding  SSSL  resonator  when  the  discontinuity  is  inserted. 
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II.  THEORETICAL  ANALYSIS 


The  microstrip  resonator  to  be  analyzed  is  shown  in  Fig.l.  A  strip 
transmission  line  which  has  a  step  in  width  is  placed  on  the  suspended  substrate. 
The  shielded  suspended  substrate  line  is  constructed  by  placing  the  substrate  within 
a  channel  milled  in  the  sidewalls.  Selection  of  both  depth  and  width  of  the  shield 
enclosure  is  important.  Dimensions  must  be  sufficiently  small  to  avoid  propagation 
of  waveguide  modes  within  the  shield:  these  modes  would  actually  produce  an 
unwanted  coupling  with  the  resonance  one  is  trying  to  detect.  On  the  other  hand, 
if  the  enclosure  dimensions  are  too  small,  the  propagation  characteristics  of  the 
transmission  line  itself  would  be  significantly  modified. 


Figure  1  End  and  top  view  of  suspended  substrate  microstrip  line 
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It  is  assumed  that  the  thickness  of  the  conducting  strip  is  negligible  and  that 
all  the  media  and  conductors  are  lossless.  For  simplicity,  the  strip  is  assumed  to 
be  located  symmetrically  in  the  longitudinal  direction  within  the  ends  of  the 
shielding  enclosure. 
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A.  FOURIER-TRANSFORMED  GALERKIN’S  METHOD  FOR  THE 


SHIELDED  STRIP  RESONATOR 

The  analysis  of  wave  propagation  on  shielded  strip  transmission  systems  of 
the  type  considered  here  has  been  carried  out  by  Itoh  and  Uwano  fRef.2].  These 
authors  treated  the  problem  using  Galerkin’s  method  of  moments  in  the  Fourier- 
Transform  domain.  Their  work  is  summarized  in  terms  of  a  set  of  Green’s 
function  equations: 


(la) 


^  =  E 


In  Eqs.  (la)  and  (lb),  and  are  the  z  and  x  components  of  the  currents  on 
the  strip  conductor,  and  and  E^  are  the  components  of  electric  field  tangent  to 
the  substrate  surface,  and  the  are  the  Green’s  functions. 

The  tildes  over  the  factors  in  Eqs.(la)  and  (Ib)  imply  Fourier  transformation 
of  the  respective  quantities.  In  this  work,  the  Fourier  transformation  is  carried  out 
in  the  bounded  region  interior  to  the  shielding  enclosure  placed  around  the  line 
segment  representing  the  resonator.  Therefore  we  have  the  finite  Fourier 
transformation: 


5 


(2a) 


^  a 

f  (ic^,  P^)  =£  ^f(x,y,  z)  e^'^'^e^^-^^dxdz 

'T  'T 


where  and  p„  are  the  discrete  transform  variables  defined  by. 


K  = 

(n-l/2)7r/a 

for  Ej  even,  odd(in  x)  modes 

(2b) 

K  = 

ntr/a 

for  Ej  odd,  -H,  even(in  x)  modes 

(2c) 

Pn,  = 

(m-l/2)7r/RL  for  E^  even,  odd(in  z)  modes 

(2d) 

P„.  = 

mre/RL 

for  E,  odd,  even(in  z)  modes 

(2e) 

where  a  is  width  of  waveguide,  and  RL  is  length  of  waveguide  enclosure. 

The  analytical  task  then  consists  of  assuming  suitable  coordinate  forms  of  the 
cunent  densities  and  J,,,  and  Fourier -transforming  these.  Inner  products  aie  then 
formed  in  Eqs  (la)  and  (lb),  in  accordance  with  the  standard  procedure  of  the 
method  of  moments.  This  procedure  gives  rise  to  a  set  of  homogeneous  equations 
in  the  unknown  coefficients  assumed  in  and  The  solution  condition  for  the 
simultaneous  equations  leads  to  a  determination  of  the  resonance  fr  equency  of  the 
strip  line  resonator  of  the  problem.  The  unknowns  and  E,^  can  be  eliminated 
by  applying  Galerkin’s  method  in  the  spectral  domain.  The  first  step  is  to  expand 
the  unknowns  and  in  terms  of  assumed  basis  functions  and  with 
unknown  coefficients  c„,  and  d„,  . 
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The  basis  functions  and  must  be  chosen  to  be  the  Fourier  transforms 
of  space-domain  functions  Jzn,(x,z)  and  3^n,{XyZ)  which  are  physically  realistic,  and 
which  are  zero  except  for  the  region  kl  <  w,  and  bl  <  f.  Substituting  (3)  into 
(1)  and  taking  inner  products  of  the  resulting  equations  with  the  basis  function  and 
for  different  values  of  i  yields  the  matrix  equations, 


(4b) 

__  w  _  ^  _ 

^  XZ^2  ^rr'^  xnT  ®  i-1 , 2  ,  .  .  ,  I'l 

m-l  ir-^l 


The  right-hand  sides  of  (4)  are  zero  by  virtue  of  Parseval’s  theorem,  because 
the  currents  JJx),  and  the  field  components  Ej(x,h2),  E,j(x,h^)  vanish  in 
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complementary  regions  of  x.  For  example,  if  the  inner  product  of  E,  on  the  right- 
hand  side  of  Eq.(l)  with  J^jCkn)  is  taken,  one  obtains: 


(x)  ( -x)  dx=0 


In  the  above,  J^,(x)  is  zero  outside  the  strip,  and  E^(x)  is  zero  on  the  strip. 
In  this  way,  the  final  boundary'  condition  is  satisfied. 

Equations  (4)  will  be  expressed  in  matrix  form  as  follows: 


N  M 


E 

i  =  l, 2,  . 

..N 

(5a) 

K 

M 

L 

i=l ; 2; . 

(5b) 

/7-: 

where  from  the  definition  of  the  inner  products  associated  with  the  Fourier 
transform  defined  by  (2),  the  matrix  elements  are:. 


K, 


,1) 


H'l  IT-l 


zz'^  zw  P  rj 


(6a) 
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(6b) 


K 


(1,2) 


n~l  m-1 


„IZzxJ  »,<■*:»•?.) 


z-(z.l) 

^2in 


n-l  /n-1 


(6c) 


n-l  ni’l 

A  homogeneous  system  of  equations  is  thus  obtained  in  terms  of  the  unknown 
coefficients  c„,  and  dn,.  In  order  that  c^,  and  dn,  have  nontrivial  solutions,  the 
determinant  of  the  matrix  must  be  zero,  and  hence  the  resonant  frequency  is 
determined  for  the  resonator  and  discontinuity  structures  assumed. 

Equations  (5)  are  now  solved  for  the  resonant  frequency  co^,  by  seeking  the 
root  of  the  resulting  characteristic  equation.  The  resonance  frequency  of  the 
suspended  stripline  resonator  is  derived  from  the  obtained  value  of  (o„[Eq.(8h)]. 

The  accuracy  of  the  solution  can  be  systematically  improved  by  increasing  the 
number  of  basis  functions  (M,N)  and  by  solving  larger  size  matrix  equations. 
However,  if  the  first  few  basis  functions  are  chosen  so  as  to  approximate  the  actual 
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unknown  current  distribution  reasonably  well,  the  necessary  size  of  the  matrix  can 
be  held  small  for  a  given  accuracy  of  the  solution,  resulting  in  numerical  efficiency. 
Hence  the  choice  of  basis  functions  is  important  from  the  numerical  point  of  view. 
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B.  GREEN  S  FUNCTION  FOR  THE  MOMENT  METHOD 

Equations  (la)  and  (lb)  may  be  expressed  in  the  matrix  form:[Ref.2] 


FJ  12..  2JF 


where  the  transformed  Green’s  Functions  have  the  values: 


sr  “  o  h'^ 


xs  ~  ^  zx~ 


^nP;r 


IZ  -Z  J 


z .  = 


_ y  ys^c.^y  yj^c, _ 

’V  c.  -‘C.  r.  -1^+a  c. 

r,  ..  o,  Y,,  Yj,- 
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Y  zz^c.*y  zi^t- 


y  ziy  z^C,C,^^y  z3C,C,^^y  ^,y  ,,C,C,^^y‘ 


(6e) 


where  the  symbols  C,.have  the  values  specified  below: 
case  >  0, 


C,^=coth  YjCf,  C,^=cothY2t,  C^^cothy^h  (8f) 


case  Y,‘  <  0. 


C,=-jcoty^d,  cot  y^t,  C^^=-jcoty,h  (8g) 


where 


Y  i=/^I+P^rM'^e7^ 


(8h) 


y 


,  Y, 

JoaCj 


{8i) 


Y  zi 


y  X 


(8j) 


where  subscripts  i  refer  to  the  i"*  medium. 

The  quantities  Z„,  and  Z„  are  actually  the  Fourier  tninsfomis  of 

dyadic  Green’s  function  components. 
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C.  CHOICE  OF  RESONATOR  FREQUENCY  BELOW  THE  CUTOFF  OF 
THE  SHIELD  ENCLOSURE 

The  cutoff  frequency  of  the  dielectric-loaded  waveguide  formed  by  the 
shielded  SSSL  shown  in  Fig.l  is  found  to  be  5.73GH2,  as  calculated  by  the 
equation  below  [Ref.5]. 


f  -  -£L  t  (gj-l)  (9a) 

^  ■  2a ^|  b  e, 

where  c  is  the  velocity  of  light  and  the  value  of  e,  is  2.35.  The  dimensions  a,b 
and  t  are  defined  in  Fig.  1,  and  numerical  values  are  shown  in  Table  I. 

TABLE  I  CONFIGURATION  DATA 


Dimension 

Data 

Dimension 

Data 

a 

22.86  mm 

t 

1.59  mm 

b 

10.16  mm 

h 

4.285mm 

5.00  mm 

Er 

2.35 

d 

1 

1 

4.285mm 

RL 

100.00  mm 

Dimensions  in  Table  I  correspond  to  Fig.l  as  follows: 
a  =  Width  of  shield 
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b  =  Height  of  shield 

d  =  Height  of  lower  air  layer 

t  -  Thickr.ess  of  dielectric  substrate  layer 

h  =  Height  of  upper  air  layer 

w,  =  Microstrip  line  width  in  continuity  case 

e,  =  Dielectric  constant  of  substrate  layer 

RL  =  Length  of  shield 

The  effective  dielectric  constant  and  the  frequency  determine  wavelength 
in  the  SSSL: 


X 


s 


(9b) 


where  c  is  the  velocity  of  light,  f  is  the  frequency,  and  e^ff  is  the  effective  dielectric 
constant.  Thus,  X^2  determines  the  fundamental  resonance  of  the  open  ended 
length  of  strip,  C. 

For  simplicity,  the  strip  is  assumed  to  be  symmetrically  located,  although  the 
present  method  of  analysis  can  be  easily  extended  to  more  general  cases.  Also, 
the  operating  frequency  is  taken  to  be  below  cutoff  of  the  shielding  waveguide 
partially  filled  with  substrate  material,  in  order  to  avoid  interaction  of  the  strip 
transmission-line  resonance  with  standing-wave  resonances  of  the  waveguide  modes 
within  the  shielding  enclosure. 
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D.  EQUIVALENT  NETWORK  FOR  A  DISCONTINUITY  IN  A 
TRANSMISSION  LINE  SYSTEM:  THE  RESONATOR  METHOD 

In  the  resonator  method  for  evaluating  an  equivalent  circuit  representation  for 
a  discontinuity  in  a  transmission-line  system,  the  discontinuity  is  placed  at  the 
approximate  center  of  an  open-ended  strip  of  the  transmission  line. 

Ignoring  the  end  effect,  the  fundamental  resonance  of  the  unperturbed 
resonator  strip  occurs  at  a  total  resonator  length  of  XJ2  of  the  guided  wave  length 
of  the  system.  For  this  mode  there  is  an  current  maximum  at  the  center  of  the 
resonator  (Fig.  2a].  The  next  higher  resonance  of  the  open-ended  resonator  occurs 
at  a  resonator  length  of  of  the  system  wavelength  [Fig.  2b]. 

In  any  event,  if  Che  line  resonator  is  housed  within  a  metallic  shield,  for 
example  of  waveguide  cross-section,  it  is  essential  to  choose  a  resonator  length 
such  that  its  resonance  frequency  is  well  below  cutoff  of  the  propagating  modes  of 
the  waveguide  system  involved.  This  is  a  dielectric-loaded  waveguide  because  of 
the  presence  of  the  dielectric  substrate. 

In  the  resonator  method  for  determining  the  equivalent  ciicuit  representation 
of  a  discontinuity,  the  discontinuity  is  assumed  to  be  lossless,  and  the  equivalent 
circuit  is  regarded  as  a  two-port  circuit  element  acting  as  an  impedance  transformer 
between  the  two  line-segments  which  make  up  the  resonator[Fig.  3]. 
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Figure  2  Resonance  modes 
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Figure  3  Impedances  for  line  segments 


Considering  the  resonator  transmission  line  segment  at  the  left  of  the  two- 
port,  the  occurrence  of  resonance  is  specified  by  the  condition: 

h  =  Zinl  (10) 

where  the  star  means  complex  conjugate,  z,  is  the  impedance  at  the  input  of  the 
open  stub  on  the  left,  and  Zj^,  is  the  impedance  seen  looking  into  the  uiput  at  port 
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1.  The  two-port  is  terminated  in  the  second  open-ended  line  segment  which 
makes  up  the  remainder  of  the  resonator[Fig.  4]. 


, _  ^inl 

L . . .  - . 

_ 1 

1  H 

— -c 

Figure  4  Equivalent-circuit  expression  for  discontinuity 


If  the  discontinuity  under  test  is  physically  symmetrical  about  its  center  point, 
then  it  may  be  represented  by  a  symmetrical  T-  or  PI-  circuit[Fig.  5).  However, 
for  an  asymmetrical  discontinuity,  such  as  for  example  the  step  ci.ange  in  width, 
the  discontinuity  may  be  represented  by  a  general  unsymmetrical  T  or  PI,  or  an  L- 
section  if  appropriate. 

In  the  evaluation  of  the  discontinuity  by  the  resonance  method,  a  pair  of  stub 
lengths  {|  and  Cj  is  chosen,  of  lengths  appropriate  to  the  test  ((|+f2)=^p/2,  or 
(C.+f;)s^g.  The  spectral  domain  model  of  the  discontinuity  is  tlien  employed  to 
calculate  the  resonance  frequency  f„  of  the  resonator  plus  discontinuity.  The 
terminating  stub  length  f,  is  changed  by  a  small  amount,  and  then  the  new  length 
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of  the  stub  at  the  input  which  is  required  to  re-establish  resonance  at  the  same 
frequency  f^  as  above,  is  determined.  This  process  is  repeated  until  a  sufficient 
number  of  stub-length  pairs  ((,,£2).  which  produce  resonance  at  the  same  frequency 
f„.  is  available  for  the  calculation  of  the  element  values  of  the  desired  equivalent 
circuit  representation. 

In  the  case  when  an  L-section  equivalent  circuit  model  is  wanted,  it  turns  out 
that,  although  only  2  independent  circuit  elements  are  included  in  the  equivalent 
circuit,  a  determination  of  3  pairs  of  (£,,£2)  values  is  required,  in  order  to  effect 
separation  of  variables  in  the  circuit  equations.  This  is  shown  in  section  III. 

Likewise,  when  an  equivalent  T-  or  Pl-section  presentation  is  evaluated  for 
an  asymmetrical  discontinuity,  a  total  of  4  of  the  (£i,£2)  pairs  must  be  resonated. 
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in  order  to  provide  the  data  needed  for  the  determination  of  the  3  independent 
circuit  elements  in  the  equivalent  circuit. 
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UI.  COMPUTER  PROGRAM  CONSTRUCTION 


A.  CURRENT  DISTRIBUTIONS 

In  actual  computations  for  the  dominant  mode,  the  strip  current  densities 
and  J„,  have  been  chosen  to  have  physically  plausible  distributions,  as  shown  in 
Fig.  6,  for  the  continuous  strip  which  has  width  of  2w,  and  length  of  2t 


1.  Continuity  Case 
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The  coordinate  forms  of  the  current  distributions  in  Fig.  6  are: 


^  2W^ 


(11a) 


(lib) 


1  TT  V 

J^(a')  =  — sin — ^ 

h\ 


(11c) 


JAz)  =  - 


(lid) 


Fourier  transforms  of  these  aie; 


2sin(A-„hi, )  3 

J ,  {k)  - _ _JLJ_  + _ I _ [XIN] 


(12a) 


/  VT7T-  ,1.  1  2sin(;c„Kj)  2  [1-cos  (ir^h',)  ) 

(  XIN  =  cos  (kw.)  ~ _ , _  4- - ^_J —  ) 

KPl 


2] 


4ncosp^i 


27tsin(/c„Vj) 


cos  (p^i) 

Pn-Z' 


sin(p^i) 

(P.i)^ 


And  then,  forming  products  of  the  x-  and  z-  dependent  factors, 


J,(A„,p;,)=c7,(;c,j-j^.(p^) 


^^.<^.^P.)=-^(ic„)-J,(p„) 


(12b) 


(12c) 


(12d) 


(13a) 


(13b) 
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2.  Step  Discontinuity  Case 


The  coordinate  forms  of  current  distributions  m  Fig.  7  are: 


2w^ 


(14a) 


fl+—  3 

2  (Vj 


(14b) 


(14c) 


T  /«\  Z 


(14d) 


J^,  ix)  =-i-sin— 

Wl  U’j 


(14e) 


J  .(x)  =  — sin  — 
u’,  tv. 


(14f) 


2  if 


(14g) 


2  if 


(14h) 


The  Fourier  transformed  distributions  of  Eqs.  05)  are; 
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^  =Same  as  Eg.  (125) 


(15a) 


(^n»^2) 


[XIN] 


(15b) 


(  ;fiN H  CO.  ^  1 1-COS  (;c„>.,)  ] 


kh\ 


(^.^z) 


'■^  2i  ( P  m^  =Saine  as  Eg.  ( I2b) 


(15c) 


47tCOS  (P„ij) 


(15d) 


Eg.  (12c) 


(15e) 


(15f) 
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(16b) 


J_,..  (A.,  p;,)  =  c,J 


ik^)  -J 


where  the  individual  terms  have  non-zero  values  only  in  the  line  segment  indicated 
by  the  corresponding  subscript,  as  explained  above. 

From  Eqs.  (5)  and  (6) 

-  -  (17a) 

E  ^  p  J  ^  zzJ  ^K’  P 

n=J  m=l 
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(17b) 


n=l  m'l 


X21‘E  E  fj',,-  p  .>  z  ..j  «  (J'r.  P  J 


(170) 


n=l  w-1 


K22^E  E'^'^i^K-PJ  2'  „  ( Jc„,  P  „) 


(17d) 


n'2  /j!=3 
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B.  COMPUTATION  PROCEDURES 


1.  Change  of  Transform  Variable 

In  the  present  work,  the  strip  line  resonator  was  fully  enclosed  in  a  metal 
shield.  This  permitted  the  use  of  a  finite  Fourier  transform  rather  than  the  integral 
transform  which  is  typical  of  the  infinite-line  calculations. 

The  inner  products  can  then  be  carried  out  as  truncated  summations, 
saving  much  computer  timefRef.9]. 

p  CB 

Jo 

Strictly  speaking,  the  limits  are  n=-<x5,  m=-o«>  on  the  low  side,  unless  symmetr>' 
holds,  and  the  sum  is: 

The  wavenumbers  used  for  the  two  summations  are: 

x-direction 

z-direction 


E 


EE 

/)*i  /p-i 
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As  shown  on  Table  II  and  III,  a  summation  over  n  of  20  terms  is  enough, 
and  for  the  summation  over  m,  100  terms  leads  to  convergence  within  1  MHz  error 
which  is  negligible  in  comparison  with  the  resonant  frequency. 


TABLE  II  VARIATION  OF  NUMBER  OF  TERMS  N 


n 

m 

"o  (GHz) 

10 

100 

5.7168 

20  * 

100 

5.7250 

40 

100 

5.7261 

TABLE  III  VARIATION  OF  NUMBER  OF  TERMS  M 


n 

m 

fo  (GHz) 

20 

40 

5.7072 

20 

80 

5.7228 

20 

100  ♦ 

5.7250 

20 

200 

5.7275 
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2.  Equations  for  Some  Factors 


a.  Propagation  Constant  P 

By  use  of  the  program  shown  in  Appendix  A,  P  is  computed  at  each 
resonant  frequency  and  line-width  w,  and  Wj.  This  program  is  constructed  using 
the  moment  method  of  calculation. 


b.  Effective  Dielectric  Constant 

The  calculation  of  uses  P(f),  frequency-dependent  propagation 
constant,  and  P„^,  quasi-static  propagation  constant  [Ref.3j. 


'  ©ff 


-  r  P  ^ 1 2 
P, 


air 


(18) 


c.  Characteristic  Impedance  Zff) 

The  computation  of  characteristic  impedance  uses  an  empirical 
calculation  for  the  quasistatic  value,  "LfO)  which  has  been  confirmed  to  be  accurate 
[Ref.  10].  This  value  is  then  corrected  to  the  wanted  frequency  f  by  use  of  e,ff  in 
the  expression  :[Ref,4] 


6f  f 


(19) 
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for  0  <  w  <  a/2; 


ZJO)  =  60[V+J^-ln[ 


1+4  ( 


V  =  -1.7866  -  0.2035(t/b)  +  0.4750(a/b) 
R  =  1.0835  +  0.1007(t/b)  -  0.09457(a/b) 


for  a/2  <  w  <  a; 

Z^(0)  =  1207t[V+i?(.^'4l.3930  +  0. 6670- ln(2l^4l. 444) 

jd  b 

V  =  -0.6301  -  0.07082(t/b)  +  0.2470(a/b) 

R  =:  1.9492  +  0.1553(tA))  -  0.5123(a/b) 
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C.  PRESENIATION  OF  THE  L-CIRCUIT  COMPONENTS  WHICH 
REPRESENT  THE  DISCONTINUITY 

The  network  is  considered  as  a  two-port,  i.e.,  4-terminal,  element  tenninated 
in  open-ended  stub  sections,  as  indicated  in  Fig.  8. 


Figure  8  L-circuit  representation 


The  impedance  at  the  input  to  port  1  is  Zj„|.  The  two  stubs  have  input 
impedances: 


=  J  ZojCO  t  P  j  i  j 

(20a) 

(20b) 

The  input  impedance  to  the  network  is  then: 
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Setting  Zjni  =  z,*,  and  expanding  the  expression,  we  get: 


(22a) 

where  we  define;  A  s  z,*  -  Z2 
Bsz,‘ 

C  s  -2,*  2, 

We  now  assume  that  we  have  values  A,  B,  C  and  A’,B’,C’,  as  well  as 
A",B".C",  of  the  constants  defined  above,  corresponding  to  different  stub-length 
pairs,  ((,  ,£2).  (£/,£20)  and  (£'i",f20,  all  of  which  have  been  determined  by 
spectral-domain  calculations  to  resonate  at  the  same  resonance  frequency  f„. 

We  rewrite  the  equation  above  for  each  of  the  stated  cases: 


z^A’-z^Z^^z^B'^d 

(22b) 

z.^A"-z^z^*z^B''-d' 

(22c) 

Eqs.(22)  are  solved  simultaneously  to  obtain  the  values  for  the  circuit 


components  z,  and  7.^: 


(23a) 


(C-C")  (B-B0-(C-C0  [B-B") 
{A-A")  -  (A-A')  {B-B") 


{C-C')-Z^kA-A') 

(B-BO 


(23b) 


The  component  values  for  a  three-element  PI-  or  T-circuit  can  be  calculated 
in  a  completely  analogous  manner,  but  with  use  of  a  set  of  four  ((,,(2)  pairs. 
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D.  RESULTS  OF  COMPUTATION 


1.  Strip  Resonator  ■with  Uniform  Width 

a.  Resonant  Frequencies  with  Varying  Height  of  Lower  Air  Layer 
The  lower  air  gap  is  varied  from  the  bottom  of  the  shield  to  the 
middle  of  the  shield  height.  When  the  dielectric  material  is  near  the  bottom, 
comparable  to  the  case  of  the  microstrip  line  configuration,  the  resonant  frequencies 
of  SSSL  show  the  same  dependence  on  width  w,  as  for  a  microstrip  resonator 
(Ref.iJ.  As  the  air  gap  becomes  larger,  the  tendency  is  reversed[Fig.9J. 
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b.  Resonant  Frequencies  with  Varying  Lengths  of  Stub  Segments  on 
the  Test  Resonator 

When  the  length  of  strip  is  greater,  the  resonant  frequency  is  smaller 
as  expected,  because  the  resonant  wavelength  is  proportional  to  the  length  of  strip. 


Figure  10  Resonant  frequency  vs.  length  of  strip  with  chajige  of  lowei  air  gap 


c.  Effects  of  Change  of  Substrate  Dielectric  Constant 

As  the  substrate  dielectric  constant  is  increased,  the  resonant 
frequency  of  the  strip  resonator  is  reduced.  But  the  variation  of  resonant 
frequency  with  change  of  lower  air  gap  is  faster  as  dielectric  constant 
increases[Fig.ll]. 


Figure  11  Resonant  frequency  vs.  width  of  strip  with  change  of 
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2.  Strip  Resonator  with  Step  Discontinuity  in  Width 


a.  Resonant  Frequencies  with  Discontinuity  at  Center  of  Resonator 
When  the  length  C|(=C2)  is  decreased,  the  resonant  frequency  is 
increased  as  in  the  uniform  width  resonator  shown  in  Fig.  10.  A  greater  difference 
between  widths  w,  and  Wj  results  in  a  larger  decrease  of  resonant  frequency. 


Figure  12  Resonant  frequency  vs.  Ci(=C2)  with  different  Wj/Wj 
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b.  Resonant  Frequencies  with  Varying  Ratio  of  Lengths 


A  change  of  the  ratio  of  the  lengths  affects  the  resonant  frequencies. 
In  turn,  the  greater  difference  between  the  lengths  of  the  segments  of  different 
widths  the  higher  is  the  resonator  frequency.  However,  the  resonant  frequency 
converged  to  a  single  value,  as  the  ratio  of  two  lengths  is  was  made  larger. 
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Figure  13  Resonant  frequency  vs.  Cj/Cj  with  different  w,/w2 
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c.  Capacitance  of  Step  with  varying  Parameters 


The  equivalent-circuit  capacitance  of  Fig,8  is  dependent  on  the 
resonant  frequency.  To  get  the  different  resonant  frequencies  the  length  of  stub, 
{i(={2)*  is  changed.  When  (2  becomes  shorter,  there  is  corresponding  C,  which 
produces  resonance  with  the  same  frequency  as  the  reference  case.  As  mentioned 
before,  3  pairs  of  ((,,  Cj)  are  needed.  Here,  the  reference  case,®, ={2,  can  not  be 
one  of  the  three  because  it  was  found  that  using  this  case  resulted  in  negative 
values  of  capacitance  resulting  from  the  computation.  More  investigation  is 
required  to  determine  the  reason  for  this  behavior. 

As  shown  in  Fig.  14  &  Fig.  15,  a  greater  value  of  the  ratio  w,/w2 
corresponds  to  a  greater  value  of  discontinuity  capacitance,  as  expected[Ref.7]. 
However,  when  the  Wj/Wj  approaches  to  unity,  the  capacitance  does  not  decrease 
continuously  but  was  found  to  increase.  Since  the  limiting  case,  w,/w2=l,  should 
correspond  to  zero  capacitance,  the  result  obtained  seems  to  be  due  to  the  fact  that 
the  difference  terms  in  Eq.(23)  represent  the  taking  of  differences  between  nearly- 
equal  quantities,  and  a  small  error  in  the  individual  quantities  will  produce  a  large 
error  in  the  result.  Other  factors  are  the  approximate  nature  of  the  basis  functions 
and/or  the  small  number  of  basis  functions. 

Shorter  values  of  f ,  and  (2  correspond  to  higher  resonant  frequency. 
The  variation  of  resonant  frequency  shows  the  frequency-dependence  of  the 
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equivalent-  circuit  capacitance  [Fig.  16].  Also,  a  greater  substrate  dielectric  constant 
corresponds  to  a  greater  equivalent-circuit  capacitance  as  expected  on  physical 
grounds  [Fig.  17]. 
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Figure  14  Capacitance  vs.  Wj/Wj  with  change  of  initial  resonator  length 


Figure  15  Capacitance  vs.  Wj/Wj  with  change  of  e, 
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d.  Inductance  of  Step  with  vaiying  Parameters 


The  procedure  for  the  calculation  of  the  equivalent  circuit  inductance 
is  analogous  to  that  used  in  the  case  of  the  c{q)acitance.  The  calculated  inductance 
did  not  converge  to  zero  when  Wj/Wj  goes  to  unity,  also[Ref.8].  The  inductance 
does  not  show  a  monotonic  incmase  for  Wj/Wj  greater  than  two[Fig.  18].  This 
seems  to  be  due  to  the  approximate  nature  of  the  calculation. 

Fig.  19  generally  shows  similar  behavior  of  the  calculated  inductance, 
with  several  different  values  of  dielectric  constant.  The  inductance  decreases  as  the 
resonant  frequency  increases,  however[Fig.20].  Increasing  the  substrate  dielectric 
constant  causes  a  uniform  increase  in  the  discontinuity  inductance[Fig.21]. 
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Figure  18  Inductance  vs.  Wj/Wj  with  change  of  initial  resonator  length 


Figure  19  Inductance  vs.  w/wj  with  change  of  e, 
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Figure  20  Inductance  vs.  resonant  frequency  with  change  of  w/wj 
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IV.  CONCLUSION 


This  thesis  follows  the  methods  originated  by  Itoh  for  the  analysis  of  the 
shielded  suspended  stripline  resonator[Ref.l  &  3].  A  fully-shielded  enclosure  was 
employed,  allowing  the  use  of  a  finite  Fourier  transform  in  two  coordinate 
directions.  This  change  reduced  the  computation  time  significantly,  while 
producing  accuracy  comparable  with  that  of  the  use  of  the  integral  transform  along 
the  line  axis. 

The  perturbed-iC'  onator  technique  permitted  the  use  of  the  strip  current- 
density  distributions  suggested  by  Itoh[Ref.l].  These  are  known  to  give  accurate 
results  in  Galerkin’s  method. 

We  may  assume  various  functional  forms  for  the  current  distribution  on  the 
line-stub  of  the  resonator.  However,  no  single  distribution  can  give  an  exact  result, 
though  a  larger  number  of  distribution  functions  should  lead  to  more  accurate 
results.  The  choice  of  distribution  function  is  important  because  of  the  very  small 
number  of  functions  used  in  this  work. 

The  data  of  Fig.9  shows  an  interesting  change  of  the  slope  of  the  resonant 
frequency  vs  width  ratio  with  varying  height  of  lower  air  gap[Fig.l].  Fig.  14  &  18 
shows  increasing  of  capacitance  and  inductance  respectively,  as  the  resonant 
frequency  decreases.  But,  as  mentioned  above,  these  graphs  show  some  deviation 
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from  expectation,  in  that  the  capacitance  and  inductance  values  fail  to  become  zero 
with  no  step  in  width.  This  behavior  should  be  given  further  research. 
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APPENDIX.  COMPUTER  PROGRAMS  USING  FORTRAN 


A.  RESONANT  FREQUENCY 


FILEt  STEPFO  FOR  A1 


INTEGER  I, N,M, COMPARE, COUNT 
REALX8  W1,W2,L1,L2,A,D,T,H,ER,F 
REALX8  KN, BETA, OMEGA, RL,RRL 

REALX8  PI,E0,MU0,C 
REALX8  GAM1S,GAM2S 

REALX8  JZIN1,J2IN2,JZKN1,JZBE1,JXKN1,JXBE1,JX1,JZ1 
REALX8  JZKN2,JZBE2, JZBE3,JXKN2, JXBE2 

REALX8  CMP, LIMIT, OLDF,NEHF,LIMF 
C 

LOGICAL  LLOG 
C 

COMPL  EXX16  GAMAl , GAMA2 , GAMAS, CTl , CT2, CT3 
C0MPLEXX16  ZEN,ZENR,ZED,ZHN,ZHD,ZHNR,ZDEN 
C0HPLEXX16  ZE, ZH, ZZZ, ZZX, ZXZ, ZXX 
C0MPLEXX16  K11,K12,K22,J,SK11,SK12,SK22,DET 
CHARACTERX8  FNAME 
CHARACTERX12  SNAME 
C 

0PEN(1,FILE='DATAIN1»,STATUS=»0LD*) 

PRINTX,  'ENTER  OUTPUT  DATA  FILE  NAME  <FNAME>=' 
READX,  FNAME 
0PEN(2,FILE=FNAME) 

SNAME=FNAME//'PLT' 

OPEN(3,FILE=SNAME) 

HRITE(2,x)  •  FILE  NAME  =  ', FNAME 
WRITE(3,X)  '  FILE  NAME  =  », SNAME 

PARAMETERS 
C=3.E8 

J=CMPLX(0. ,1 . ) 

PI=4.XATAN(1 . ) 

E0  =  1 .E-9/(36  .XPI) 

MU0=4.XPIX1 .E-7 
C  LIMIT=l.E-2 

LIMIT=l.E-5 
C0MPARE=1 
C0UNT=1 
LL0G=.TRUE. 

C 

CALL  READF(W1,W2,L1,L2,A,D,T,H,ER) 

PRINTX, 'ENTER  THE  NUMBER  OF  OUTER  SUM  LOOP  N  =  » 
READX, N 
C  N=20 

WRITE(2,x)  'THE  NUMBER  OF  OUTER  SUM  LOOP  N  =',N 

WRITE(3,X)  'THE  NUMBER  OF  OUTER  SUM  LOOP  N  =»,N 

C 

PRINTX,  'ENTER  THE  NUMBER  OF  INNER  SUM  LOOP  M  =  ' 
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oo  o  o  o  o  r>  o  o  t-o 


READX,  M 

C  M=20,  40,  100,  500,  1000,  2000,  3000,  4000  ETC. 

WRITE(2,X)  'THE  NUMBER  OF  INNER  SUM  LOOP  M  =',M 

WRITE(3,X)  'THE  NUMBER  OF  INNER  SUM  LOOP  M  =‘,M 

C 

RL=0.10D0 

RRL=RLX1.D3 

C 

HRITE(2,x)  *RL  =  SRRL 
WRITE{3,X)  »RL  s  SRRL 
C 

1  PRINTx, 'ENTER  INITIAL  FREQUENCY  (GHZ)  • 

READX,F 
F=F*1.E9 

PRINT*, 'FREQUENCY  =  ',F 
WRITE(2,X)  'FREQUENCY  =  •,F 

0MEGA=2.*PI*F 

K11=CMPLX(0.,0.) 

K12=CMPLX(0.,0.) 

K22=CMPLX<0.,0.) 

DO  20  1  =  1, N 
KN=I*PI/A 
KN=(I-.5)*PI/A 

JZIN1=C0S(KNKW1)-2.XSIN(KNKW1)/(KNXW1) 

+  +2.X(1-C0S(KNXW1))/(KNXW1)XX2 

JZIN2=C0S(KNXW2)-2.XS1N(KNXW2)/(KNXW2)+2.X(1-C0S(KNXN2)) 
+  /(KNXW2)XX2 

JZKN1=2.XSIN<KNXH1)/(KNXW1)+3.XJZIN1/(KNXH1)XX2 
JZKN2=2.XSIN(KN*W2)/(KNXW2)+3.XJZIN2/(KNXW2)XX2 
JXKN1=-2.XPIXSIN(KNXW1)/(PIXX2-(KNXW1)XK2) 
JXKN2=-2.XPIXSIN(KNXW2)/(PIX*2-CKNXW2)XX2) 

SK11=CMPLX<0.,0.) 

SK12=CMPLX(0.,0.) 

SK22=CMPLX(0.,0.) 

DO  40  IM=1,M 

BETA=(REAL(IM)-.5)XPI/RL 
BETA=REAL(IM)XPI/RL 

GAM1S=KNXx2+DETAxx2-EOxMU0XOMEGAxx2 

CALL  GAMCT(GAMA1,GAM1S,CT1,D) 

CALL  GAMCT(GAMA3,GAM1S,CT3,H) 

GAM2S=KNXX2+BETAXX2-EO*MUOXERXOMEGAXX2 

CALL  GAMCTCGAMA2,GAM2S,CT2,T) 

ZEN=GAMA2XCT1/ER  +  GAMA1XCT2 
ZENR=-J/(OMEGAXE0) 

ZED=CT1XCT2+CT1*CT3XGAMA2/(GAMA3XER) 

+  +CT2XCT3  +  ERXGAMA1/GAMA2 

ZE=ZENRXZEN/'ZED 
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oo  oo  DO  o(M  on 


ZHN=GAMA1XCT1  +GAMA2XCT2 
ZHNR=JXOMEGAXMUO 

ZHD=GAMA1XCT1KGAMA2XCT2  +  GAMA1KCT1KGAMA3XCT3 
+  +  GAMA2XCT2XGAMA3XCT3  +  GAMA2XX2 

ZH=ZHNRXZHN/ZHD 
C 

ZDEN=KNXX2  +  BETAXX2 
ZZZ=-(ZExBETAXK2  +  ZHXKNXX2)/ZDEN 
ZZX=-KNXBETAX(ZE-ZH)/ZDEN 
ZXZ=ZZX 

ZXX=-(ZEXKNXX2  +  ZHXBETAXX2)/ZDEN 
C 

JZBE1=A.xPIxC0S(BETAxL1)/(PIXx2-(2.xBETAxL1)xx2) 

JZBE2=A.XPIXC0S(BETAXL2)/(PIXX2-(2.XBETAXL2)XX2) 

J21=JZKN1XJ2BE1+JZKN2XJZBE2 

JZ1=JZ1X0.5D0 

JXBEl=C0S<BETAxLl/2.)/{BETAxLl/2. ) 

+  -  SIN(BETAXLl/2.)/(BETAxLl/2.)xx2 

JXBE2=C0S(BETAXL2/2. )/(BETAXL2/2. ) 

+  -  SIN<BETAXL2/2.)/(BETAxL2/2.)xx2 

JX1=JXKN1XJXBE1+JXKN2XJXBE2 

C 

SK11=SK11+JZ1XZZZXJZ1 
SK12=SK12+JZ1XZZXXJX1 
SK22=SK22+JX1XZXXXJX1 
AO  CONTINUE 

WRITE(2,X)  JZ1,JX1 

K11=K11+SK11 
K12=K12+SK12 
K22=K22+SK22 
CONTINUE 

PRINTX,  'CHECK  Kll  =  ',K11 
WRITE(2,X)  'CHECK  Kll  =  •,K11 
PRINTX,  'CHECK  K12  =  ',K12 
WRITE(2,X1  'CHECK  K12  =  ’,K12 
PRINTX,  'CHECK  K22  =  ',K22 
WRITE(2,X)  'CHECK  K22  =  ’,K22 


DET=K11XK22-K12XK12 
PRINTX,  'DET  =  '/DET 
HRITE(2,xi  'DETERMINENT  =  ',DET 


CALL  FDS0RT(F, DET, COUNT) 
C0UNT=C0UNT+1 


IF  (LLOG)  THEN 

PRINTX,  'MAKE  SURE  OPPOSITE 
PRINTX,  'IF  SAME  SIGNS  THEN 


SIGNS  BETNEEN  1ST  8  2ND  DETS' 
ENTER  9  ,  OR  ENTER  O' 
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oo  o  oo 


READX,  UK 
IF  (IJK.EQ.9)  THEN 
0LDF=F 
GOTO  1 
ENDIF 

LL0G=. FALSE. 

ENDIF 


CMP=REAL(DET) 

LIMF»ABS(F-OLDF) 

PRINTX,  'CHECK  LIMIT  FREQ  »,LIMF 
WRITE(2,K)  'CHECK  LIMIT  FREQ  =  '^IMF 

IF  (ABS(CMP).GT. LIMIT. AND. LIMF.GT.1.E3)  THEN 
IF  (ABS(CMP).GT. LIMIT. AND. LIMF.GT.l.E-1)  THEN 
CALL  SLVAL(CMP,OLDF,F/NEWF, COMPARE) 

OLDF=F 
F=NEWF 
GOTO  11 
ENDIF 

PRINT*,  'LAST  FREQUENCY  =  ',F 
WRITE(2,K)  'LAST  FREQUENCY  =  ',F 

CALL  FDSORT<F,DET,100) 

STOP 

END 


SUBROUTINE  READF(W1,W2,U,L2,A,D,T,H,ER) 
REALX8  Wl,W2,Ll,L2,A,D,T,H,ER 
INTEGER  DDDD 

PRINTx,'DATA  INPUT  i  KEYBORD(I), FILE(O) ' 
READX,DDDD 

9  IF(DDDD.EQ.O)  THEN 

READ(1,X)  W1,W2,L1,L2,A,D,T,H,ER 
ELSEIF(DDDD.EQ.l)  THEN 

PRINT*, 'ENTER  THE  HI  8  W2' 

READ*, HI, H2 

PRINT*, 'ENTER  LI  &  L2' 

READ*, LI, L2 
PRINT*, 'ENTER  A' 

READ*, A 

PRINT*, 'ENTER  D,  T,  8  H» 

READ*,D,T,H 

PRINT*, 'ENTER  ER' 

READ*,ER 

ELSE 

PRINT*, 'SELECT  I  OR  0  ' 

READ*,  DDDD 
GO  TO  9 
ENDIF 

PRINT  32,  W1,W2,L1,L2,A,D,T,H,ER 
WRITE(2,32)  W1,W2,L1,L2,A,D,T,H,ER 
WRITE(3,32)  W1,W2,L1,L2,A,D,T,H,ER 
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W1  »  WlKl.E-3/2. 

H2  =  W2X1.E-3/2. 

LI  =  LlXl.E-3 
L2  =  L2X1.E-3 
A  =  AKl.E-3/2. 

D  =  Dxl.E-3 
T  =  Txl.E-3 
H  =  HKl.E-3 

32  F0RMAT(T3, 'UNITtMILIMETER' 

+  ,/,T3,'Wl  =  •,F9.5,5X, •M2  =  •,F9.5 

+  ,/,T3,'Ll  =  •,F9.5,5X,*L2  =  »,F9.5 

+  ,/,T3,'A  =  •,F9.5 

+  ,/,T3,'D  =  •,F9.5,5X,'T  =  •,F9.5,5X,'H  =  ',F9.5 

+  ,/,T3,'ER  =  •,F9.5,/) 

RETURN 

END 

CXXXXXKXXXKXXKXXXXiOtXKXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXX 

SUBROUTINE  GAMCTCGAMA, GAMS, CT, LENGTH) 

C0MPLEXX16  GAMA,CT 

REALX8  GAMS, LENGTH, GAMR,CTR 

GAMR  =DSQRT(DABS(GAMS)) 

IF(GAMS.LT.O.DO)  THEN 

GAMA  =  CMPLX(0.D0,GAMR) 

CTR  =  -l./TAN(GAMRXLENGTH) 

CT  =  CMPLX(0.D0,CTR) 

ELSE 

GAMA  =  CMPLX(GAMR,0.D0) 

CTR  =  l./TANHCGAMRXLENGTH) 

CT  =  CHPLX(CTR,0.D0) 

ENDIF 

RETURN 

END 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

SUBROUTINE  SLVAL ( FUN, OXV, XV, NEWXV, SN) 

INTEGER  SN 

REALX8  FUN,OXV,XV,NEMXV,POSXV,NEGXV 

IF(SH.EQ.l)  THEN 
IFCFUN.GT.O. )  THEN 
POSXV  =  XV 
NEGXV  =  OXV 
ELSE 

NEGXV  =  XV 
POSXV  =  OXV 
ENDIF 

NEHXV  =  (POSXV  +  NEGXV)/2. 

SN  =  SN  +2 
ELSE 

IF(FUN.GT,0. )  THEN 
POSXV  =  XV 
ELSE 

NEGXV  =  XV 
ENDIF 

NEMXV  =  (POSXV  +NEGXV)/2. 

ENDIF 

RETURN 

END 


CXXXXXXXXXXXXXXXXXXXXXXXXXKXXXKXXXXXKXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXX' 

SUBROUTINE  FDSORK  FQY, SDET, ICOUNT) 

INTEGER  ICOUNT/IC 

REALXB  FREQY(100),FQY,RTEMP,MAGDET,TFQY 
C0MPLEXX16  DETMT(100),SDET,CTEMP,TSDET 
TFQY  =  FQY 
TSDET  =  SDET 
IF(ICOUNT.EQ.l)  THEN 
FREQY(l)  =  TFQY 
DETMT(l)  =  TSDET 
ELSEIF(ICOUNT.NE.IOO)  THEN 
DO  10  I=1,1C0UNT-1 

RTEMP  =  FREQY(I) 

CTEMP  =  DETMT(I) 

IF  ( RTEMP. GT. TFQY)  THEN 
FREQY(I)  =  TFQY 
TFQY  =  RTEMP 
DETMTd)  =  TSDET 
TSDET  =  CTEMP 
ENDIF 

10  CONriNUE 

FREQY(ICOUNT)  =  TFQY 
DETMTdCOUNT)  =  TSDET 
IC  =  ICOUNT 

ELSEIF(ICOUNT.EQ.IOO)  THEN 
WRITE(3,110) 

WRITE(3,120) 

no  F0RMAT(//,T8, 'FREQUENCY  *,T23,'MAG  OF  DET  ' , TA3, 'DETERMINENT' ) 

120  FORMAT<1X,TAO,'RE. •,T56,'IM.*) 

DO  20  J  =1,IC 

MAGDET  *DSQRT(REAL(DETMT<J))XX2+DIMAG(DETMT(J))XX2) 

HRI TE<  3,210)  FREQY( J ) , MAGDET, DETMT( J ) 

210  FORMAT<T5,F15.2,T20,F15.3,T38,2(E15.8)) 

20  CONTINUE 

ENDIF 
RETURN 
END 
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B.  PROPAGATION  CONSTANT  P 


FILE.  BETA  FOR  A1 


INTEGER  MMAX,NMAX 

REALX8  N.A,D,T,H,PI/ER/F6,OMG,EO,KN 

REALMS  G1S2,B0,B00,DBTA,FLAG»AA2,CC,WN 

REALXB  G1,CT01,CT03,GT1,GT3,G2S2,G2,CT02,KE1,KH1,J21 

REALX8  B,OUT,EREFF,FC,BIM,Q,R 

DIMENSION  B(60),OUT(30) 

C0MPLEXX16  GM1,GM3,CT1,CT3,GM2,CT2,DENE,ZE,KE,DENH,ZH,KH 
COMPLEX  TRMN,SUMUP/JZ11 
C 

0PEN(2,FILE=*BETAVAL*,STATUS=‘0LD») 

C 

PI=4.D0XATAN(1.D0) 

FLAG=0 

C 

MMAX=30 

NMAX=30 

C 

PRINTX, ‘ENTER  FREQ., GHZ.' 

READJ<,  FG 
C 

PRINTX, ‘ENTER  SHLD.  WDTH(MM)' 

REAO)(,  AA2 
A=AA2/2000 

PRINTX, 'ENTER  SHLD.  HT.(MM)‘ 

READX,CC 
CC=CC/1 . D3 

PRINT*, ‘ENTER  LOWER  AIRGAP,(MM)* 

READ*,D 

D=D/1.D3 

PRINT*, 'ENTER  SUBSTRATE  THICKNESS(MM) ' 

READ*,T 

T=T/1.D3 

H=CC-(D+T) 

PRINT*, 'ENTER  THE  SUBSTRATE  DIELEC.  CONST.  ' 

READX,ER 

PRINT*, 'ENTER  THE  LINEWIDTH(MM) » 

READ*,WW 

W=WW/2.D3 

Q=0 . 2077+1 . 2177XT/CC-0 . 0836AKAX2 ./CC 
PRINT*, '0=',Q 

R=0. 03451-0. 1031*T/CC+0.01742*A*2./CC 

PRINT*, 'R=',R 

TT=Q-RKL0G(W*2./CC) 

SS=LOG(l./SQRT(ER)) 

EREFF=1./(1+TT*SS) 

EREFF=EREFF**2 
PRINT*,  'EREFF  =  ',EREFF 


0MG=FGXPI*2E9 
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E0=8.8<iD-12 

C  BO«(OMG/3D8)KSQRT((ER+l)/2) 

BO=(OMG/3D8)KSQRT(EREFF) 

BOO-O.2OD0 
DBTA^O.IODO 
GO  TO  4 
5  BOBBIN 

BOO=0.85DO 
DBTA’^O.OOTDO 
A  DO  2  M=1,MMAX 
SUMUP=0 

B(M)>BOX(BOO-t-MxDBTA) 

DO  1  N=1,NMAX 
KN=CN-0.5D0)KPI/A 
G1S2=KNXX2+B(M)XX2-0MGX)(2/9D16 
G1=SQRT(ABS(G1S2)) 

IF(GlS2.LT.O)  THEN 
GM1=CMPLX(0.D0>G1) 

GM3SCM1 

CT01*-1.DO/TAN(01XH) 

CT03=-1.D0/TAN(G1XD) 

CT1=CMPLX(0.D0,CT01) 

CT3=CMPLX<0.D0,CT03) 

ELSE 

GM1=CMPLX(G1,0.D0) 

GT1=1.D0/TANH(G1XH) 

GM3=GM1 

GT3=1.D0/TANH(G1XD) 

CTl=CMPLX<GTl,O.DO) 

CT3=CMPLX(GT3,0.D0) 

ENDIF 

C 

G2S2=KNXX2+B<M)XX2-ERX0MGXX2/9 . D16 
G2=SQRT(ABS(G2S2)) 

IF  (G2S2.lt. 0)  THEN 
GM2=CMPLX(0.D0,G2) 

CT02=-1.D0/TAN<G2XT) 

CT2=CMPLX(0.D0,CT02) 

ELSE 

GM2=CMPLX(G2,0.D0) 

GT2=1.D0/TANH(G2XT) 

CT2=CMPLX(GT2,0.) 

ENDIF 

C 

DENE=CT2XCT3+CT1XCT3X(GM2/GM1)/ER+CT1XCT2+{GM3/GM2)XER 

ZE=(GM2XCT3/ER+GM3XCT2)/DENE 

KE1=-1 .DO/(OMGXEO) 

KE=CMPLXCO.DO,KEI) 

ZE=KEXZE 

DENH=GM1XCTIXGM2XCT2+GM1XCTIXGM3XCT3+GM2XCT2XGM3XCT3+GM2XX2 

ZH= ( GM2XCT2+GM3XCT3 )/DENH 

KH1=0MGXPIXAD-7 

KH=CMPLX(O.DO,KHI) 

ZH=ZHXKH 

JZ1=COSCKNXW)-2.DOXSIN(KNXW)/(KNXW)+2.DOX(1.-COSCKNXW))/(KNXW)XX2 

JZI=2.D0XSIN(KNXH)/(KNXW)+3.D0/((KNXW)XX2)XJZ1 

JZII=CMPLX<JZ1,O.DO) 

TRMN=-(B(M)XX2XZE+KNXX2XZH)/(B(M)XX2+KNXX2)J;JZ11XX2 
SUMUP=TRMN+SUMUP 
1  CONTINUE 
C 
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OUT<M)sCABS(SUMUP) 

CONTINUE 


DO  3  M=2,(MMAX-1) 

IF  <(OUT(M).LT.OUT(M-l)).AND.(OUT(M).LT.OUT(M+l)))  THEN 
BIN=B(M) 

ELSE 

ENDIF 

CONTINUE 

IF  CFLAG.EQ.O)  THEN 
FLAG=1.D0 
GO  TO  5 
ELSE 

PRINTX, 'SUSP.  STRIPLINE  BETA  =  SBIN 
WRITE(2,X)  ‘SUSP.  STRIPLINE  BETA  =  •,BIN 
ENDIF 
STOP 
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C.  EQUIVALENT  CIRCUIT  CAPACITANCE  &  INDUCTANCE 


FILE.  EQICIR  FOR  A1 


INTEGER  I 

REALX8  BETA1,BETA2»Z01,Z02 
REALX8  PI, F, OMEGA, Ml, M2 

COMPL  EXX16  A ,  B ,  C,  J ,  Z3, 24,  NUMZ5,  ;,ENZ3,  NUMZ4 ,  DENZ4 
C0MPLEXX16  CAPAC,INDUC,MC,ML 
DIMENSION  A(3),B(3),C(3) 

CHARACTERX8  FNAME 

PRINTX,  'ENTER  OUTPUT  DATA  FILE  NAME' 

READX, FNAME 
0PEN(2,FILE=FNAME) 

J=CMPLX(0.D0,1.D0) 

PI=4.D0XATAN(1.D0) 

PRINTX, 'ENTER  RES.  FREQUENCY  (GHZ)' 

READX,F 

NRITE(2,X)  'RESONANT  FREQUENCY(GH2)  =',F 

F=FX1.E9 

OMEGA=2.DOXPIXF 

PRINTX, 'ENTER  Ml  8  M2' 

REA0X,M1,M2 
MRITE(2,130)  Ml, M2 

130  F0RHAT(T3,'M1  =  ' , F3 . 1, 'MM* ,5X, 'M2  =  ',F6.4,'MM') 
PRINTX,  'ENTER  ZOl  8  Z02' 

READX,  201,202 
PRINTX, 'ZOl  =  •,Z01 
PRINTX, *Z02  =  ',Z02 
MRITE<2,100)  201, Z02 

100  FORHAT(T3,'Z01  =  ' , F8 . 4, 5X, »Z02  =  »,F8.4) 

PRINTX,  'ENTER  BETAl  8  BETA2’ 

READX,  BETAl, BETA2 
PRINTX,  'BETAl  =  ', BETAl 
PRINTX,  'BETA2  =  ',BETA2 
MRITE(2,110)  BETAl, BETA2 

no  F0RMAT(T3, 'BETAl  =  ',F8 .4,5X, 'BETA2  =  ',F8.4) 

DO  10  1=1, 3,1 

CALL  Z1Z2ABC(BETAI,BETA2,Z01,Z02,J,I,A,B,C) 

10  CONTINUE 

NUMZ5=(C(1)-C(3))X(BC1)-B(2))-(C(1)-C(2))XCB(1)-B(3)) 
DENZ3=<A(1)-A(3))X(B(1)-B(2))-(A(1)-A(2))X(B(1)-B(3)) 
NUMZ3=(C(1)-C(2))X(B(2)-B(3))-(C(2)-C(5))X(B(1)-B(2)) 
DENZ3=(A(1)-A(2))X(B(2)-B(3))-CA(2)-A(3))X(B(1)-B(2)) 
Z3=NUMZ3/DENZ3 

NUMZ4=(C(1)-C(2))-Z3X(A(1)-A(E)) 

DENZ4=B(1)-B(2) 

Z4=NUMZ4/DENZ4 

PRINTX, 'NUMZ3  = ' , NUMZ3, 'DENZ5  =',DENZ3 
PRINTX, »NUMZ4  =',NUMZ4,*DENZ4  =»,DENZ4 
MRITE(2,X)  'NUMZ3  =',NUMZ3, 'DENZ3  =',DENZ3 
MRITE(2,X)  'NUMZ4  =',NUMZ4, 'DENZ4  =',DENZ4 
PRINTX, »Z3  =  ',23 
PRINTX, *Z4  =  ',24 
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WRITE(2,X)  »23  =  ',23 
WRITE(2,K)  'Z^i  =  ',24 
CAPAC=l./fZ3XJX0ME6A) 
INDUC=ZA/(JXOMEGA) 

PRINTX,  'CAPAC.s  •,CAPAC 
PRINTX,  'INDUC.s  •,INDUC 
WRITE(2,X)  'CAPAC.  =  •,CAPAC 
WRITE(2,X)  'INDUC.  =  *,INDUC 
WC=1.DO/(OMEGAXCAPAC) 
NL=OMEGAXINDUC 
PRINTX,  'l/WC  =  •,WC 
PRINTX,  'WL  =  ',WL 
WRITE(2,X)  »1/WC  =  •,HC 
WRITE(2,X)  'WL  =  ',WL 
STOP 
END 


SUBROUTINE  Z1Z2ABC(BETA1,BET A2,Z01,Z02, J, I , A, B, C) 
INTEGER  I 

REALX8  BETA1,DETA2,Z01,Z02,L1,L2 
COMPLEXXl^  A,B,C,J,Z1,Z2 
DIMENSION  A<3),B<3),C(3) 


PRINTX,  'ENTER  LI  8  L2  (MM)' 

READX,  L1,L2 
PRINTX,  'LKMM)  =  ',L1 
PRINTX,  »L2(MM)  =  ',L2 
HRITE(2,120)  L1,L2 

120  F0RMAT(T3, 'LKMM)  =  » ,  F8 . 4, 5X, » L2(MM)  =  ',F8.4) 
L1=L1/1.D3 
L2=L2/1.D3 

Z1=JX201/TAN<BETA1XL1) 

Z2=-JX202/TAN(BETA2XL2) 

PRINTX, »21  =  ',21 
PRINTX, '22  =  ',22 
WRITE(2,X)'21  =  ',21 
WRITE(2,X)  '22  =  ',22 
A(I)=C0NJG(21)-22 
B(I)=C0NJG<Z1) 

C(I)=-C0NJG(Z1)X  22 

PRINTX, 'A( ',1, '  )  =  »,A(I) 

PRINTX, 'B( ',1, '  )  =  »,B(I) 

PRINTX, 'C< ',1, '  )  =  ',C(I) 

HRITE(2,X)  »A(',I,»  )  =  »,A(I) 

WRITE(2,X)  'BC,!,'  )  =  ’,B(I) 

HRITE(2,X)  'C(',I,'  )  =  *,C(I) 

RETURN 

END 
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